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Abstract 

Cellular responses often require the fast activation or repression of specific genes, which depends 
on Transcription Factors (TFs) that have to quickly find the promoters of these genes within a 
large genome. Transcription Factors (TFs) search for their DNA promoter target by alternating 
between bulk diffusion and sliding along the DNA, a mechanism known as facilitated diffusion. 
We study a facilitated diffusion framework with switching between three search modes: a bulk 
mode and two sliding modes triggered by conformational changes between two protein confor¬ 
mations. In one conformation (search mode) the TF interacts unspecifically with the DNA 
backbone resulting in fast sliding. In the other conformation (recognition mode) it interacts 
specifically and strongly with DNA base pairs leading to slow displacement. From the bulk, a 
TF associates with the DNA at a random position that is correlated with the previous disso¬ 
ciation point, which implicitly is a function of the DNA structure. The target affinity depends 
on the conformation. We derive exact expressions for the mean first passage time (MFPT) to 
bind to the promoter and the conditional probability to bind before detaching when arriving at 
the promoter site. We systematically explore the parameter space and compare various search 
scenarios. We compare our results with experimental data for the dimeric Lac repressor search 
in E.Coli bacteria. We find that a coiled DNA conformation is absolutely necessary for a fast 
MFPT. With frequent spontaneous conformational changes, a fast search time is achieved even 
when a TF becomes immobilized in the recognition state due to the specific bindings. We find a 
MFPT compatible with experimental data in presence of a specific TF-DNA interaction energy 
that has a Gaussian distribution with a large variance. 

Keywords: Facilitated diffusion; Transcription factor; Mean first passage time; Gene regula¬ 
tion; Mathematical model; Lac repressor; E.Goli 
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Introduction 


Transcription factors (TFs) regulate gene activation by binding to DNA promoter sites. To 
enable a fast cellular response that relies on the activation or repression of specific genes, TFs 
perform a facilitated diffusion search where they alternate between three dimensional (3D) dif¬ 
fusion in the bulk and ID diffusion (sliding) along the DNA (for reviews see ([ll-@)). Initially, 
facilitated diffusion was introduced to explain the experimental finding that the in-vitro associ¬ 
ation rate of the Lac-I repressor with its promoter sites placed on A-phage DNA was around 100 
times larger than the Smoluchowski limit ~ for a 3D diffusion process (0). Theo¬ 

retical considerations showed that a search that alternates between 3D diffusion and ID sliding 
can have a higher association rate compared to a pure 3D search QM)- With dilute DNA the 
search time is dominated by the 3D excursions between subsequent DNA binding events, and 
sliding increases the association rate by enlarging the effective target size (antenna effect). Later 
on single molecnle techniques provided a direct experimental proof of the facilitated diffusion 
mechanism (UMi). 

It has also soon been realized &M) that frequent bindings to the DNA are problematic 
becanse sliding along the DNA is slow due to strong TF-DNA interactions ([IJ, Il4l ). In a 
dense DNA environment with frequent bindings to the DNA and slow ID diffusion, the antenna 
effect becomes negligible and facilitated diffnsion is slower compared to a pure 3D search. For 
example, in E.Coli with a volume IV_ ~ the measured search time of the Lac repressor for 

its promoter site is r ~ 350s 0, , 'X'liis corrGsporicis to ciii OiSSOCiSitioii rOito 

10®M“^s“^, much lower than the Smoluchowski limit. If a TF could specifically bind only to its 
promoter site and bounce off from the rest of the DNA, the search time would be extremely fast 
around ~ V/{^iiDR) ~ 5s. However, because a TF cannot already recognize its target from the 
bulk, frequent DNA associations are essential. Thus, the question arises: How is a fast search 
possible within a large genome despite of facilitated diffusion ? 

When a TF is bound to the DNA and interacts with the underlying base pairs (bps), the 
diffusion coefficient for sliding decays exponentially with the variance of the binding energy 
distribution SS). With a simple facilitated diffusion model that comprises sliding along the 
DNA and uniform redistributions in 3D, one finds that a search time of the order of minutes is 
only compatible with a variance < l.^ksT (fl^ . In contrast, bind ing energy estimates for the 
Cro and PurR TF reveal a much larger variance around 5 — tSksT (1^, 1^). This indicates that 
a simple diffusion process is not sufficient to explain the search dynamics when a TF is bonnd 
to the DNA. It has been proposed that a TF switches between two protein conformations with 
different binding affinities to the DNA SEl). In the search conformation, a TF interacts only 
non-specifically with the DNA backbone leading to a smooth energy profile and fast diffusion. 
In the recognition conformation, a TF interacts specifically with the underlying DNA sequence 
resulting in a rough energy landscape and slow diffusion. Conformational changes of the TF 
protein are indeed supported by experimental observations (12, 2fll- 241. 

In this work we investigate a general framework for a facilitated diffusion search with con¬ 
formational changes. We analytically derive the mean first passage time (MFPT) to bind to 
the target and the conditional probability to bind before dissociation when a TF arrives at the 
target site. We further compute the ratio of the time spent in the bulk compared to attached to 
the DNA, the apparent diffusion constant for sliding along the DNA, and the average sliding dis¬ 
tance before detaching. We consider a search process with Poissonian switchings between three 
states (Fig. [T|). State 1 and 2 (recognition and search mode) correspond to two different protein 
conformations with conformation dependent TF-DNA interactions. Therefore the diffusion co- 
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Figure 1: Facilitated diffusion framework. (A) Schematic of the search process: In state 3 
a TF is freely diffusing in the bulk. It attaches to the DNA at a random position following a 
Gaussian distribution centered around the previous detaching position. In state 1 and 2 a TF 
is attached to the DNA and diffuses along the DNA. In state 1 it specifically interacts with the 
underlying DNA sequence and diffusion is slow. In state 2 it non-specifically interacts with the 
DNA backbone and diffusion is fast. The binding affinity to the target depends on the state. 
(B) Switchings between states occurs with Poissonian rates. The rate ki 2 depends on the energy 
of the specific interaction. 


efficients for sliding along the DNA and the target affinity depend on the conformation. In state 
3, a TF is diffusing in the bulk and it associates with the DNA at a random position following a 
Gaussian distribution centered around the previous dissociation point. By modifying the target 
affinity in state 2 we evaluate the impact of induced switchings at the target site and the effect 
that a TF misses the target when arriving at the target site. By varying the correlation distance 
between dissociation and association point we estimate how the DNA conformation and coiling 
affect the MFPT. With our analytic expressions we can precisely evaluate the whole parameter 
space. We analyze various search scenarios within the same framework, which is important 
to accurately compare results. Other approaches partly rely on MFPT analysis, kinetic the¬ 
ory, thermodynamic equilibrium considerations, scaling arguments and other approximations, 
which complicates a comparison of results obtained with different methods and approximations 
0, l25l-[^). A clean MFPT analysis for a 3 states switchin g ni odel with maximal target affinity 
in state 1 and no affinity in state 2 has been performed in (31). Gompared to (01. the authors 
additionally consider 3D excursions using a closed-cell approach as described in (j^). However, 
due to the difficulty to compute the 3D kernel, only the asymptotic limits corresponding to 
uniform redistributions and a rod-like DNA are discussed. The first passage time distribution 
for a switching process between two ID states has been studied in (j33lb 


Model 


Model description and MFPT analysis 

We start by presenting the mathematical framework and the MFPT analysis. We postpone the 
description of the biological motivation to the results part. We consider a search that switches 
between 3 states with Poissonian switching rates kij (Fig. [1]). In state 1 and 2 a TF is attached 
to the DNA of length 2L and slides with state dependent diffusion constants Di and ZI 2 . To 
simplify the analysis we consider that the target is located at the center. An off-centered target 
results in a higher MFPT up to maximally a factor of 4 if the target is located at the periphery 
(assuming that the MFPT scales ~ (3J, l35l ) . When a TF reaches the target it binds with 

state dependent affinities xi and X 2 - In state 1 a TF switches to state 2 with rate ki 2 - In state 
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2, in addition to switching to state 1 with rate k 2 i, a TF can also dissociate with rate ^23 and 
switch to state 3 where it diffuses in the bulk. From state 3 it associates with the DNA with 
the rate ^32 at a random position drawn from a Gaussian distribution with variance cj| centered 
around the previous dissociation point. 

Because of the Gaussian attaching distribution, we model state 3 as an ID diffusion process 
along the DNA with an effective diffusion constant D 3 = ^4^ and no target affinity (x 3 = 0). 
Thus, we finally arrive at a framework with switchings between three ID states. The backward 
Fokker-Planck equation for the probability p{x,t,n\y,m) to find the TF at time t in state n at 
position X, conditioned that it started at t = 0 in state m at position y, is (j36l . 13711 


dtpix,t,n\y,m) = Dmdyp{x,t,n\y,m)-2xmP{x,t,n\y,m)5{y) 

3 

- ^ kmi {p{x, t, n\y, m) - p{x, t, n\y, i)) , 
i=l 

with reflecting boundary conditions at y = ±L. The mean sojourn time spent in state n is 

roo pL 

Xn,m{y) = / dt dxp{x,t,n\y,m). 

Jo J-L 

From Eq. [T] we find that the Tn,m{y) satisfy the system of equations 


( 1 ) 


( 2 ) 


dJm’Tji fjiiy) 'y ^ ikm+dmi kmi) Tn^i{y) ‘JXmd{y)Tn^m{y) — 


(3) 


i=l 


with km+ = In the Supplementary Information (SI) we exactly solve Eq. [3] and 

derive analytic expressions for the Tn^miy)- The MEPT when initially in state m at position 
y is Tm{y) = Yln=i'''n,m{y)■ We focus on the MEPT with uniform initial distribution, fm = 
^ j_i^Tm{y)dy. Because switchings between states occur fast compared to the overall search 
time, the dependency of fm on m is negligible. Furthermore, the mean sojourn times fi^m 
approximately satisfy the scaling relations (see SI Eq. 40) fi^m '■ T 2 ,m ■ Tz,m 
Hence, the MFPT with uniform initial distribution is well approximated by 


_ . ^21 . ksi 

fcl2 ■ fc23 ' 




^21 K21K32/ 


1 1 /C23 1 

kl2 k2l k2l ^32 


(4) 


where N 12 = Ti,ifei 2 is the average number of switchings between states 1 and 2 . 

To reveal the scaling of the MEPT as a function of the DNA length L, we introduce the 
scaled length L = ^ defined with the reference length Lq = Ibp. To simplify the discussion, we 
focus on a scenario with maximal affinity in state 1 , xi = 00 , in which case case the search is 
over when the TF encounters the target in state 1 (the analysis in the SI is performed with a 
general xi). For a long DNA {L'^pi S> 1) we compute in the SI 


Ni2 = LI 


12 


( 


Cl,l 


Cl,2 


\yiVJLl P2y/P2j 


+ L 


2 ^ 12 ^ 21^32 


3,5 


(5) 
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with the following dimensionless parameters: 

r2^.. 

riQivij 


kj — jj , ci,i — hih 2 ( ——I- bi 


aid 


I 0‘2d 

Cl,2 = 121^32 ( — -1- 02 


a = 


^12 + ^21 + hs + ^32 j /3 — ^21^32 + ^12(^23 + ^32) , hl/2 ~ 2 ^ -\/a^ — 4 


hl{h2 + hi — h2) h2ih2 + hi — hi) 

oi — — i — r -^—N—5 ®2 — 


^ 12(^1 “ h2) 


K2 = 


L 0 X 2 

D 2 


d = 


^ 12(^2 — hi) 
1 


hi = 


h2 — hi 
h2{hi - M 2 ) 


h2 = 


h2 — h2 
h2{h2 - hi) 


h2K2 


hi 62 


y/hn y/h2 


ai 02 /21 

e = H-— + 


yfhl y/h2 112^2 


Search with uniform redistribution in state 3 and no target affinity in state 2 

With uniform redistributions (cjs = oo) we have ^32 = 2Lq/(t| = 0. With ^2 = 0 and ^32 = 0, 
Eq. [5] simplifies to 




hi - h2 y/hl h2 - hi y/Th 


( 6 ) 


in agreement with (j32l ). We note that the term ~ in Eq. [5] vanished and we now have W 2 ~ E, 
which leads to a faster search for large L. As stated before, Eq. [ 6 ]is valid for L^hi (* = 1)2) 
and therefore cannot be applied for ^23 ^ 0. For /c 23 —^ 0 one eigenvalue, say /ri, vanishes and 
the condition L^hi S> 1 is violated. For example, with a fixed L and /C 23 0 the TF remains 

bound to the DNA. In this case we expect that the MEET scales ~ L? and not ~ L, which is 
indeed the case, as can be shown by a refined analysis. However, for any fixed value ^23 > 0, by 
increasing L, N 12 eventually scales ~ L due to the uniform redistributions. 


Optimal switching scenario in state 2 

When the properties of state 1 and state 3 are fixed (and D 2 is fixed), we compute the optimal 
switching rates /c 2 i and ^23 that minimize the MEET. We introduce the parameters 


af = 


2Di 

ki2 


cj| = 


2 D 2 


k21 + ^23 


q = 


k23 


k21 + k23 


c = 


d. 

< 7 ^ 


(7) 


af is the mean square displacement in state 1 , ctI is the mean square displacement in state 2 
before switching either to state 1 or 2 , is the detaching probability, and ( is the ratio of the 
displacements in state 1 and 2. We use the variables q and ( instead of k 2 i and ^ 23 . We have 
^23 = h2QC and /21 = ^ 2(1 — q)C- Because diffusion in state 1 is slow compared to state 2 , we 
have C ^ 1- We further consider that the probability to switch from state 2 to 1 is much larger 
than the dissociation probability, such that g <C 1. For C ^ 1 and g <C 1 we obtain from Eq. 0] 
and Eq. [ 6 ] the asymptotic 



err 


ki2 2E)2C 


+ 


k32 


( 8 ) 


For fixed ai, D 2 , ki 2 and /C 32 , the minimum of r with respect to (g, () is 

y/2L 


(Tlkl2 


(1 + 7C)^ 


r = 


(9) 
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achieved at 7C = and 'y'^qC = 1, where 7 = and 6 = ^7 = Interestingly, 

whereas the optimal rate k 2 i depends on the properties of state 1 , we find for /c 23 the optimal 
value ^23 = ^qC, = "§ 7 ^ ~ ^32) independent of the properties of state 1 . 

Search with two states only 

To derive the MFPT with switchings between two states we set ^23 —)• 0. With ki = 00 and 
^2 = 0 (state 2 now corresponds to the bulk state without binding) we find 


with = h 2 + hi- For ki 2 ^ 0 we recover r = With uniform redistributions in state 2 
we get (/21 = 0 ) 


1 1 

ki2 k2i 


( 10 ) 


- — ft ( hi L ^ li2 1 
T — LI 12 — H — 
\h2^ h2y/h2, 


= V2 


L 

(Tl 


1 1 

ki2 k2i 


( 11 ) 


Eq. [TT] as a function of ki 2 has a minimum for ki 2 = k 2 i (compare with ^23 = k ^2 obtained with 
3 states). 


Results 

We present results that we compare to experimental measurements for a dimeric Lac repressor 
search in E.Coli bacteria. We use a DNA length 2L = 4.8 x 10® bps (jl3), a TF attaches to the 
DNA from the bulk after an average time k ^2 — l-4:ms (fl^ . 38), and the diffusion constant in 
state 2 is L )2 = 2fj,m'^fs {D 2 = 1.7 x lO’^^) (flU . We keep these values fixed throughout 
the following analysis and we focus on investigating the impact of the remaining parameters. 


Search scenario with conformational changes 

A TF is freely diffusing in the bulk (state 3) and attaches to the DNA with a Poissonian rate 
k'i 2 (FigHJ. We consider that the association position follows a Gaussian distribution with 
variance <73 centered around the previous dissociation point. <73 is the correlation distance 
between subsequent detaching and attaching positions. Hence, 73 is an effective parameters 
that implicitly depends on the DNA configuration and on coiling. For example, a uniform 
re-attaching distribution obtained for <73 = 00 is usually attributed to a highly packed DNA 
conformation (1^, ^,[^, 3^41). By varying 173 we can investigate how the DNA conformation 
affects the MFPT. We assume that a TF switches between a stable and an unstable protein 
conformation. The lifetime of the unstable conformation is short such that a TF quickly 
returns to its stable conformation after a spontaneous conformation change. When a TF is 
attached to the DNA and in the stable conformation (state 2) it non-specifically interacts with 
the DNA backbone and diffuses in a smooth potential well with non-specific energy E^s and 
fast diffusion constant D 2 . In state 2 a TF can either dissociate from the DNA with rate /C 23 , or 
switch to the unstable conformation (state 1) with rate k 2 i- The unstable conformation allows 
for additional specific TF-DNA interactions that modify the residence time k^,^ in state 1. We 
use the Arrhenius like relation A:i 2 = where AE = Ens — E (in units of ksT). We use 
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a maximal target affinity in state 1, xi = oo, in which case the search is finished when a TF 
reaches the target in state 1 for the first time. In state 2, the outcome at the target depends on 
the affinity X 2 (respectively the dimensionless parameter K 2 ). The search is over for K 2 = 00 , 
in which case a TF has maximal target affinity already in state 2. In the opposite case ^2 = 0 
a TF has no indication in state 2 that it has reached the target site and there is a probability 
that he misses the target and detaches without binding. Because switching to state 1 at the 
target site ends the search, by varying K 2 we can explore the impact of induced switchings at 
the target site. 

We use constant switching rates k 2 i, ^23 and ^ 32 . This is valid for a spatially homogenous 
DNA and a homogenous non-specific interaction in state 2. In contrast, A:i 2 and ui depend on 
the specific binding energy E and are therefore not constant along the DNA. To account for 
this, we first derive results using a constant E corresponding to a homogenous DNA, and in 
a subsequent step we average using a Gaussian distribution for E. Because of strong specific 
interactions, we focus on displacements ai that are small. The lower bound for cJi is reached 
when a TF becomes immobilized in state 1. However, this does not correspond to ui = 0, 
because a TF at least scans the base pair it binds to. A non-zero ci also accounts for stochastic 
fluctuations in the DNA position due to the switching process. By noting that the MSD of the 
maximum displacement of a diffusion process is 2 (T^, we use '^1 = to model the limiting case 
where a TF becomes immobilized in state 1 and scans only a single base pair. 

To facilitate the comparison with experimental data, we introduce the following parameters 
that characterize various properties of a search process: 

_ 1 , ^21 1 ^ _U ^ .^2 _ _2 1 I _21 “ 'I' n _ ^‘dna n 

^dna — 7 “T 7 7 ? ^Id2td — ^32^dna i ^dna — ^2 ' ’ ^dna — r> * 

k 23 k 23 kl2 q q 2 Tdna 

Tdna is the average time a TF stays bound to the DNA before detaching; rid 3 d is the ratio of 
the time bound to the DNA to diffusing in the cytoplasm; is the mean square displacement 
along the DNA before detaching; Ddna is the effective diffusion constant for sliding. 


Search with uniform redistributions and no target affinity in state 2 


We start by analyzing search processes as a function of the switching rate A:i 2 with no target 
affinity in state 2 (k 2 = 0) and uniform redistributions in state 3 (Fig. [2]). We write ki 2 as 
function of the binding strength, ki 2 = and plot quantities as a function of AE. We 

use the basal rate ^ = 10 ^s“^, which is similar to the attempt frequency used in (jlil b or 

from (39|). At this stage the exact value of ^ is not important to show the behaviour as a 
function of AE. For example, a smaller value for ^ would shift the origin of the AE-axis to the 
right in Fig. [2] and Fig. [Sj but otherwise does not affect the graphs. Later on we will estimate a 
more appropriate value for ^ by considering a Gaussian binding energy distribution. 

We compare optimal and non-optimal searches for = "^ (^ ^p is scanned in state 1) and 

ai = -^ (up to 3 bps are scanned in state 1). The optimal search is characterized by k 23 = k ^2 
and a rate k 2 i that is a function of AE and ui (see Eq. [9]). For non-optimal searches we use 
^23 = ^32 and a rate ^21 that is independent of AE, since the properties of state 1 should not 
affect the switching rate A: 2 i in state 2. We use the optimal value for A: 2 i computed with AE = 5 
(hence, k 2 i = 1.72 x for cji = ^ and ^21 = 9.92 x for ui = ^), which gives 

a fast search also with large AE. We use two different values for ai = -^ and to 

facilitate the comparison between optimal and non-optimal curves: in this case the optimal and 
non-optimal MFPT coincide for AE = 5 (Fig. [5]A). 
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Figure 2: Search with no target affinity in state 2 ( k 2 = 0) and uniform redistributions 
(fja = oo). (A) MFPT for optimal and non-optimal searches as a function of the specific binding 
strength for two values of ai {ki 2 = ^ For ci = il^/2 a TF is immobile 

in state 1 and scans only a single bp, for ui = 3/\/2 it scans three bps. The optimal MFPT 
is computed with Eq. [6l the non-optimal with Eq. [5j The blue line corresponds to the MEPT 
for a two-states switching process between state 2 and 3 with an absorbing target in state 2 
(k 2 = oo)- Energies are in units of ksT. The rest of the parameters are: ki = oo, k ^2 — l-4ms, 
D 2 = 2fim?/s, k 23 = k^ 2 - The value of A: 2 i for the non-optimal search equals the optimal value 
computed with A^; = 5. (B) Ratio of the time spent associated with the DNA compared to 
freely diffusing in the bulk. (C) Apparent diffusion constant for sliding along the DNA. (D) 
Probability P to bind to the target before dissociation when arriving at the target site. 


Fig. [2]A shows that even an immobilized TF in state 1 can have a MFPT that is compatible 
with the experimental finding 350s sm). The MFPT is faster for ai = -^ because more 
DNA is scanned during the same residence time in state 1. Interestingly, the MFPT varies only 
very little as a function of AE up to values AE ~ 5 (Fig. [2]A), suggesting that the search is 
insensitive to a large part of binding energy fluctuations. The asymptotic value of the optimal 
MFPT for small AE corresponds to a two states process where a TF switches between state 2 
and 3 and has maximal target affinity in state 2 (Fig. [2]A, blue curve). This is consistent with 
results from (j37l ) showing that a switching process can have a fast MFPT even if the searcher 
can only bind in the slow state. 

For an optimal search process with only two states (bulk and one sliding state), a TF spends 
an equal amount of time in the bulk and associated to the DNA. This is not any more the case 
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for an optimal three states process (Fig. [2jB). Although we have ^23 = /C 32 , which is similar 

to the condition for a two states process, because of switchings to state 1, a TF spends much 

more time bound to the DNA. For an optimal search we have = 1 + V^, where 6 is 

defined after Eql9j For AF = 6 and cji = -^ we obtain ri^sd ~ 10, which is similar to in vivo 

findings that a dimeric Lac repressor spends 90% of the search time bound to the DNA (ILlI ). 

The effective diffusion constant for sliding D^na decreases as specific binding becomes stronger 

(Fig. [2p). For an optimal search we compute Ddna ~ D 2 /'f'id 3 d- Thus, for ridsd = 10 we obtain 
2 

Ddna ~ 0.2^F^, which is around 4 times larger than values estimated from single molecule 
tracking experiments on flow stretched DNA (|l3l L However, such a value is in good agreement 
with results from molecular dynamics simulations for a Lac dimer (j^ L and with an apparent 
ID diffusion constant Dg// ~ 0.4^^ estimated in 0 ). Moreover, a large range of variability 
is observed for the ID diffusion constant of a Lac repressor on elongated DNA estimated from 
single molecule imaging techniques (Il4li . Whereas Ddna strongly depends on ki 2 , the sliding 
distance ddna is not affected by the residence time in state 1 (if ui remains unchanged) and 
is determined by diffusion in state 2 and the detaching rate ^ 23 - For an optimal search with 
^23 = ^ 32 , by neglecting ui, we obtain ddna ~ ~ 2206p. This value is much larger than in 

but compatible with a value around 


vivo measurements for a Lac dimer around 40 bps (I 111 . 

240 bps obtained from molecular dynamics simulations (j42lL Similar to Ddna, also for ddna a 
large experimental variability is observed using single molecule imaging techniques (Il4l ). 

Finally, when arriving at the target site in state 2, a TF can as well detach without binding 
to the target (fill, liil). To characterize such events, we compute the conditional probability P to 
bind before detaching when arriving at the target site in state 2 (see SI Eq. 61). The probability 
depends on the affinity K 2 , for example, for K 2 = 00 we have P = 1. For ^2 = 0 the probability 
is not zero and it depends on the local switching dynamics. In general, P can be expressed as a 
function of the sliding distances independent of the switching rates. Thus, for constant d\, P is 
independent of A:i 2 or AE (Fig. [2p). For the optimal search process P depends on AE because 
A: 2 i and therefore d 2 vary with AE. 


Search with finite redistributions and induced switchings at the target site 

We proceed and study the effect of the DNA configuration and induced switchings at the target 
site by varying CJ 3 and K 2 - We consider a non-optimal search with ui = ^. Lowering the 

correlation distance CJ 3 up to CJ 3 ~ ^ (around 1 % of the genome is correlated) has only little 
impact on the MFPT (Fig.[3]A.,D), but a further decrease strongly increases the MFPT (Fig. [3p). 
In contrast, induced switchings at the target site (k 2 > 0) only moderately reduce the MFPT 
(Fig. [3]A.,B). For example, the difference in the MFPT between K 2 = 0 and ac 2 = 00 is much 
smaller compared to CJ 3 ~ and CJ 3 ~ ^ (Fig. [3fl). We conclude that a reduced coiling cannot 
be compensated by induced switchings at the target site. Although a larger K 2 increases the 
probability to bind to the target (Fig. EP), this has only a minor impact because P has already 
a value around 50% for K 2 = 0. 


Search with a Gaussian binding energy distribution 


So far we used a constant rate A;i 2 corresponding to a constant binding energy E. In reality, 
A:i 2 depends on the DNA sequence and therefore on the DNA position x. To account for this, 
we consider a search with a Gaussian binding energy distribution p{E) (0, 19). We further 
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Figure 3: Search with K 2 > 0 and us < oo. We consider the search from Fig. [2]with ai = l/\/2 
and modify K 2 and (A) MFPT for various K 2 and a^. (B) MFPT with cjs = L/300 and 
various K 2 - (C) Probability to bind to the target before dissociation when arriving a the target 
site for various K 2 (note that P is independent of a^). (D) MFPT as a function of us for K 2 = 0 
and AE = 5. 


consider the case where a TF is immobile in state 1 such that ui = l/\/2 and the number of 
switchings ^"12 (see Eq. ED are both independent of x. Let w{x) be the weight function that 
measures how often position x is visited during a search process compared to the average. For 
a uniform initial distribution, from symmetry considerations, we can deduce that w{x) = 1. In 
this case the MFPT is 


T = Ni2 


w{x) 1 fc 23 1 

kl2{x) ^ k2l k21 k32 


= iVl2 


I 1 fc23 1 

ki2iE) k2i k2i /C32 


Hence, we can compute the MFPT with the average switching rate 


= 


p{E) 


ki2{E) 




E — En 


-dE. 


(13) 


(14) 


With Ens ~ —11 and a Gaussian distribution p{E) with variance <7 = 5 (Il8l . Il9l ) we obtain 


^12 ~ with C = ^ + Ens ~ 1-5. Next we checked how the results in Fig. [2] with AE = 
C = 1.5 comply with the experimental data. We hnd that a much better level of agreement 
is achieved for AE ~ 5.5. We note that the range of the AE axis in Fig. [2] depends on the 
value of Had we used a different value ~ 10^s“^, the origin of the 

AE axis would be shifted to the right and the results for AE = 5.5 in Fig. [2] would appear at 
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AE = Thus, by using ^ ~ 10 ®s“^ instead of ^ ~ 10 ^s“^ we obtain a search scenario that 

is compatible with experimental data even in presence of a Gaussian energy distribution with 

large variance. Moreover, for ^ ~ 10 ^s“^ we also have ^ ~ ^21 and we have the Arrhenius like 

relation ^ « ^-(E-Ens) ^ 
k2i 

Discussion 

We investigated a framework for facilitated diffusion with switchings between three states: a 
bulk state (state 3 ) and two states with sliding along the DNA (state 1 and 2 ) motivated by 
two TF protein conformations. The TF-DNA interaction and the target affinity depend on 
the conformation. From the bulk, a TF associates to the DNA with a Poissonian rate and a 
Gaussian distribution centered around his previous dissociation point. We analytically computed 
the MFPT, and the conditional probability to bind to the target before detaching when arriving 
at the target site. We further defined and computed various other properties that characterize 
the search process, e.g. sliding length, effective ID diffusion constant or ratio of the time spent 
in ID compared to 3 D. We compared our results with experimental data for the dimeric Lac 
repressor search in E.Goli bacteria. We investigated various properties of a search process that 
we now discuss in more detail. 


Impact of the DNA conformation 

It is still largely unclear how strongly the DNA conformation affects the search time (H-Q. In 
the literature one can find analytic results for two opposite cases: a rod-like DNA or a maximally 
coiled DNA where the re-attaching distribution is uniform. However, a systematic and consistent 
analysis where the impact of coiling is gradually changed is still outstanding. In our model, the 
association rate /C32 and the correlation distance are two effective parameters that implicitly 
depend on the DNA conformation. For fixed <73, a higher attaching k2,2 decreases the MFPT, but 
only up to a lower limit that is attained for instantaneous jumps ^32 = 00 (Eq. 0 ]). The MEPT 
is minimal for a uniform redistribution (1T3 ^ 00, Pig. [ 3 ] 4 ), a scenario that is freque ntly used to 
analyze a facilitated diffusion process with a highly packed DNA conformation ( 18 . [ 32 . 35 , & 
4 l|). We find that around 1 % of the DNA has to become correlated by the 3 D excursions in 
order to maintain such a fast MFPT (Fig. [ 3 ]A). At lower correlation distances the search time 
is greatly prolonged (Fig. [ 3 p). The value of also determines how the search time scales as 
a function of the DNA length L. To show this we consider the number of switchings N12 that 
are necessary to find the target. A12 increases proportional to L for <73 —>■ 00 (Eq. [ 5 ]), and such 
a linear dependency is usually assumed in the literature. However, for finite <73, the leading 
order asymptotic for large L is ^12 ~ and not N12 ~ L. Por finite L, a careful analysis is 
needed to determine whether the contribution ~ L or is dominant. Por our analysis we 

considered that <73 and ^32 are independent parameters, however, in general their values will be 
correlated. Por example, lets consider stretched DNA. In this case, by assuming a correlation 
distance <73 ~ L, we find N12 ~ L. However, because the DNA is stretched and a TF is diffusing 
with diffusion constant ZI3 in the bulk, we must at least have k^2 ~ ^3 1 ^ 3 - Finally, this leads to 
a MFPT that scales ~ and not ~ L. On the other hand, with strong coiling one might have 
as ^ L with a fast rate ^32 that is almost independent of <73, such that the MFPT scales ~ L. 
We conclude that without coiling it is not possible to have a MFPT that scales ~ L. Coiling is 
permissive to obtain at the same time a large correlation distance 173 and a fast attaching rate 
k 32 - 






Facilitated diffusion with three search modes 


12 


Impact of induced switchings 


If switchings from state 2 to state 1 are induced at the target site, the spontaneous switching 
can be reduced leading to a faster search because a TF spends more time in the fast state 
2 ([l^, l26l b However, it is unclear which physical mechanism would provide such a specificity. 


We investigated the impact of induced switchings by varying the target affinity in state 2 (^ 2 ). 
We find that a MFPT compatible with experimental data can be achieved without induced 
switchings (Fig. ^2]A. and Fig. We estimated that in this case a switching rate around 

k2i ^ is necessary. This implies conformational changes in the submillisecond range, that 

have also been suggested in (jl2l ). The rate ^21 could be further reduced by assuming a larger 
sliding distance cJi. If spontaneous switchings to state 1 are fast, the conditional probability P 
to bind to the target before detaching is already large for K 2 = 0 (Fig. ^SlC). In such a case 
additional induced switchings (k 2 > 0) do not much affect the search time (Fig. [3j3). Clearly, the 
impact of induced switchings would be much larger if P would be small for K 2 = 0. For example, 
this could be achieved by lowering the target affinity in state I (ki < 00 ), or by reducing the 
switching rate A: 2 i. With such conditions the MFPT would be strongly increased by blocking 
induced switchings. In general, for large K 2 = 00 the fastest MFPT is achieved by simply not 
switching to the slow state I (A; 2 i = 0). But in this case we return to a two states model with a 
bulk and a fast sliding state, incompatible with a large binding energy variability. 


Search in presence of a Ganssian binding energy profile 


State I is characterized by the displacement ai and the residence time . We analyzed the 
limiting case where a TF becomes immobilized in state I such that it scans only the base pair 
to which it binds to (ui = I/\/2). The rate A:i 2 = depends on the specific energy E 

in state I and the non-specific energy Ens in state 2. With ^ Ens = —llksT and a 

Gaussian distribution for E with variance a = bksT we obtained a MFPT around 5-6 minutes, 
compatible with in vivo experimental data for the dimeric Lac repressor IH). It is found that 
the Lac repressor dimer stays bound to the promoter for an average time around 5 minutes 
(1471) . In our model this would correspond to a target energy E = Ens — ln(CT;,) ~ —28kBT, 
compatible with data (E,[ii). At strong noncognate DNA sites with E ~ —25kBT, a TF would 
be trapped only for a short time k ^2 — ~ I2s, which resolves the trapping problem 

([ 3 ^ . The speed-stability paradox strongly relies on the assumption that a TF is found with 
high probability bound to the target at thermodynamic equilibrium. However, when the MFPT 
is of the order of minutes, such a high probability implies that a TF blocks the promoter for a 
very long time. This would impede a fast cellular response, and generate the opposite problem 
of how a promoter can get rid of a tightly bound TF. 


Conclusion and prospects 

In this work we presented a MFPT analysis for a facilitated diffusion search process with switch¬ 
ings between three states: a bulk state and two sliding states where the TF is attached to the 
DNA. The model is microscopically motivated and describes the local dynamics using effective 
parameters. Parameter values have to be extracted from more detailed models of the TF-DNA 
interaction, or by fitting our analytic expressions to experimental data for the dimeric Lac re¬ 
pressor. We focused on a qualitative analysis of the model and we showed that the model 
predictions account for many features that are observed experimentally. 
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A major simplification of the current model is the fact that we reduce the impact of the 
3D dynamics to a Poissonian association rate /C 32 and a Gaussian re-attaching distribution with 
width (T 3 . However, this simplification allowed to derive analytic results, which are important 
to precisely analyze the parameter space. Furthermore, we generalized results with uniform 
redistributions corresponding to (T 3 = 00 . Assuming that CJ 3 is correlated to the amount of DNA 
coiling, we could systematically investigate the impact of the DNA conformation. However, 
polymer models show that the distribution of CJ 3 is not a Gaussian but decays like a power law 
at large distances HHlii). The re-entry distribution is also more complicated than a single 
exponential. It will be interesting to investigate in future work how more accurate assumptions 
for the 3D dynamics based on polymer models change the results presented here. Another 
interesting project is to compute the MFPT with Levy flights in state 3 HSi). 

Instead of using a single state for the 3D dynamics with complex distributions for attaching 
time and position, one could break down the 3D dynamics into several states with simpler 
distributions and enlarge the current model by additional states. Each state would account for 
different properties of the search process, e.g. hoppings, jumps, intersegment and intersegmental 
transfers. 
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Supplementary Information 
Derivation of sojourn times and MFPT 

We start from the equations for the sojourn times 

3 

^mTn,m,iy) ^ ^ ‘^Xm^{y yo)Tn,m{fj) — ^nm (1^) 

i=l 

with the switching matrix {km+ = ^mj) 

/ h2 

Kmi — km+^mi kmi — I ^21 

V 0 

and reflecting boundary conditions at y = ±L. Because the target is located in the center at 
j/o = 0 we can restrict the analysis to the region 0 < y < L. By integrating Eq. [15] around y = 0 
we obtain (rn,m(y) = Tn,m(-2/)) 

^rnTn^rn{y)\y=0+ ~ XmTn,m{0) (17) 

which are partially reflecting boundary conditions. We remove the killing term in Eq. [15] and 
replace it with these partially reflecting boundary conditions. We introduce the dimensionless 
position X = ^, the diffusion rates Um = the dimensionless parameters Km = the 

scaled switching rates Imi = and the scaled switching matrix Lmi = Kmi/i^m- The 

scaled sojourn times 

satisfy the system of equations equations (0 < x < 1) 

3 

^ f-'mi T'n,i{x) — ^nm (t^) 

i=l 

with reflecting conditions r^^(l) = 0 at x = 1, and partially reflecting conditions f!^mi^) — 
at X = 0. In state 3 we have a reflecting boundary condition at x = 0 and x = 1. 
The functions {fn,m = Jq fn,m(x)dx) 

Vn,m{f^') — fn,m (^t)) 

have zero mean and satisfy the system of equations 

~ ^rni Vn,i{x) = —V„ m(0) ■ ( 21 ) 


— ^12 0 \ 

A:21 + k23 —k23 I (16) 

— k 32 k 32 ) 


The matrix Lmi is singular and one eigenvalue is zero. The left eigenvector to the zero eigenvalue 
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From Eq. [19] we obtain 


3 3 

m=l m=l 

and after integration we find 

3 

^ ^ fmVn,m{x) = fndi^) ~ 9n{x) i (23) 

m=l 


where 


g{x) 




(24) 


With Eq. |23| we express Vn, 2 (x} — Vn, 3 {x) as a function of Vn,i{x) and Vn, 2 {x) and then obtain 
closed system of equations for Vn,iix) and Vn, 2 ix)- By introducing 


Wnix) = Vn,l{x) - Vn, 2 {x) 


(25) 


we find from Eq. [21] 

l{x)y _ ^ hn,l{x)\ ^ ^ <,l(0) \ 

\Wn{x) J \Wn{x J \w'^{0) + gn{x)J 

with 



and 


a — li 2 + I 21 + I 23 + ^32 ; /3 


hih2 + ^ 12(^23 + ^ 32 ) • 


(26) 


(27) 


(28) 


We solve these equations as a function of ^ (0) and w'^ (0) and then compute v'^ ^ (0) and w'^ (0) 
using the boundary conditions. The eigenvalues and eigenvectors M of are 


1^1/2 = ^ (a± -4^) , 

with /3 = 9 i 92 and a = fUi + ^ 2 - With the expansion 

^ u (^x)ei + U2(x)e2 

^ Wn[X J 


(29) 


(30) 


and 


13 ^ 13 M ^ 1 -1 F2 ^ 

\0j 112(92 — 91) ^ 112(91 — 92) W / 91 — 92^ 92 — 9 i 

we derive from Eq. [20 


u'Kx) - 9 iui(x) 



/3 

112(92 — 91 ) 


+ (Wn(^) + 9n(x)) 


91 


91 - 92 


(31) 
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The equation for U 2 {x) is obtained by interchanging /xi and fj, 2 - The function ui{x) = ui{x) — 

9n{x) 

Ml—M2 


satisfies u'Kx) — /iiui(x) = —Cn,i with 


/3 


With ui{x)dx = Jq ui(x)dx = 0 we find 


hi _|_ fn 


hi - h2 hi- h2 


(32) 


(I cosh(^l-x))\ gr,{x) 

Ulix) = Cn,l - ^ + 

\hi ^Jh{smh^f\I{ ) 


hi - h2 


(33) 


The solution U 2 {x) is obtained from Eq. [33]by interchanging {cn,i,hi) with {cn, 2 ,h 2 ), where Cn ,2 
is defined in Eq. [32] with and fj ,2 interchanged. From Eq. [SOlwe obtain 


fi 2 / 1 cosh(^(l-x))\ fi 2 / 1 cosh(^(l-x)) 

Vn,l{^) = Cn,l -- ^ , -^^|+Cn ,2 - (- 


hi \hi y/hfsinh 
h2fn 

/3 


Wn{x) = Cn,l( — - 

\hl 


-g{x) 

1 cosh(y/7T(l “ x)) 


y/JIf sinh 

From Eq. [53] we further find (with fn = -j^) 


h2 \h2 y/^Z^sinhy/^ 


1 cosh(^(l - x)) 
Tt 2 y/JI^ sinh y/JI^ 


(34) 

(35) 


/ \ / \ /3 , s ^32 / N fndix) 

Vn,l{^) - Vn,3{x) = r-;—Vn,l{x) - —Wn{x) H- - - 

tl2t23 ‘23 ‘23 

Vn, 2 ix) - Vn,3ix) = Vn,l{x) - Vn,3{x) “ Wnix) 


(36) 


We complete the analysis by computing the values of v'^ ^^(0) and v'^ 2 ( 0 )- With = 0 

and ( 5 ^( 0 ) = -fn) 


~ fn 


m=l 


we can express v'^ 2 (0) as a function of v'^ ^ (0) 


<2(0) = -|^<i( 0 ) + 

‘12 ‘32 


From Eq. [32] and Eq. [38] we get 


with 


Cn,i = v'j^i{0)ai+fnbi 

Cn,2 = v'n^i{t))a2 + fnh2 


hi{h2 + hi — h2) , h2 — hi 

“1 — 1 —7-7 ) ai — 


^ 12 ( 1^1 — h 2 ) h 2 {hl — h 2 ) 

h2{h2 + hi — hi) , h2 — h2 
, t>2 = 


h2ih2 — hi) 


h2{h2 - hi) 


(37) 


(38) 


(39) 


02 = 


( 40 ) 
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To derive an equation for ^^(0) we compute rc„(0) using Eg. [351 and the relation 


'^71,1 '^71,2 — 

obtained by integrating Eq. [TUI We find 


h2 h2 


/ri\ Cn I Cn 2 * * W,l(^) /i 

W^n(O) = ^ H-^ - Cn,2^2 = -^ 

hi h2 h2 P 


- W ~ — Cn,2?2 


- N <l(0) <, 2 ( 0 ) , <l(0) 

^^^n(0) = Vi(0) - Tn,2(0) - {jn,! - Tn,2) = -1-;- 

Kl K2 112 112 


where we used 


and introduced 


From Eq. [42] we find 


Cn,l ^ Cn ,2 
hi h 2 


<,l(0) fn 

h2 fd 


coth ^/hl coth y/h^ 

si j - ? ^2 


y/hl 


y/Fl 


Cnoil + Cn,2?2 = “' 


n,l(0) '^n,2(0) fn ^n,l 


Ki 


+ 


K2 


P ^ ll 2 


1 _ _|_ 
^ 32«2 Id) 1x2 


and with Eq. [38] we obtain 

Cn,lf,l + Cn,2?2 = + /n ^ 

By inserting c„,,i and Cn,2 from Eq. [39] we obtain 

^n,l(0)(®l?l + 02^2) + fn{blf,l + &2?2) = ~^n,l(0) ( - f 1 -^ + /n f 

\Ki I 12 K 2 J V 

From this we finally get 

— fn- + - 

’ e li 2 e 


( 41 ) 


(42) 


(43) 


(44) 


(45) 


(46) 


1 1 

h2K2 fd) h2 


n,l 


(47) 


with 


d = -^ - (&16 + ^26) 

t32«^2 P ^ ^ 

e = ai 6 + 026 H-h 7 -^ . 

Kl I 12 K 2 


(48) 


We get the final expressions 
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Sojourn times and MFPT 

The scaled sojourn times are 


+ fn,l = - U„,i(0) + fn,l(0) = Un,l(x) - Un,l(0) + 


<l(0) 


Ki 


fn,2{x) = Vn,lix) - (Un,l(x) - Vn,2{x)) + fn,2 = fn,l{x) - Wn{x) + fn,2 “ fn,l 

^ , , , , ^ (50) 

= Tn,l{x) - Wn[x) ^ -;- - - 

n2 tl2 


fn,d,{x) = fn,2{x) - (VnM^) “ ^ 

t32 

where we used fn ,3 — fn ,2 = 7 ^ obtained by integrating Eq. [T9j The sojourn times with uniform 


initial distributions are 


I32 


Tn,! = -Wn,l(0) H- 

Ki 

^^n,l(0) '^( 1 , 1 ( 0 ) 


"^71,1 — 


+ 


5 '^71,2 — 'T~n,l 


<l(0) Sn,l 

h2 lu 

<,l(0) Sn,l 


- - . L/ri.,1 X X . ^n,6 /ri\ 

Tn,2 = Tn,l H-;- : - , Tn,3 = Tn,2 + 7 - (51) 

t32 


' ^nh2 ^12 

The MFPT with uniform initial distribution in state m is 


Xn,3 — Xn^2 T 


^11,3 
I 

Sn,3 


h 


32 


(52) 


r(m) = ^ t'n ^fn,m = ^ Tn,m ■ 

n=l n=l 


(53) 


Because switching between states is fast compared to the overall search time, the mean sojourn 
times are almost independent on the initial state, fw i ~ t„ 2 ~ 3 - By noting that 

we find 


T, 


t,m 


T- 




fj 


(54) 


The expression for the MFPT simplihes to 


I - I - - A I ^12 ki2k23\ 

Tl,l + r2,l + T3 ,i = Tip 1 + --h -—-— 

V A:2i k2ik32j 


= Ni2 


1 1 A:23 1 \ 

- -1 — -^23 

kl2 k21 k2lk32) 


— — + — + —^ 
k23 kl2 k23 k32) 


(55) 


where we introduced the mean number of switchings between state 1 and 2 resp. 2 and 3 

f^23 


A^12 = ^1,1^12 = fxal 


12 


iV- 


12 


k23 

k2i 


(56) 


A ^12 and N 23 can be expressed as a function of the DNA length L and the mean square displace¬ 
ments (Ti, CJ 2 and (T 3 . 
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Dependence of the search time on the DNA length 

To analyze how the MFPT depends on the DNA length L, we introduce L = where Lq = Ibp 

" L^k- ■ 

is a reference length, and extract the dependency on L. We therefore replace lij = ^ with 

where lij = is now evaluated with Lq. We proceed similarly with all the other 

parameters: Vi —)■ Ki Lk^, pi —)• a —>■ L'^a, /3 —>■ T^/3, fn o^i Oi and 

bi —>• L~‘^bi. In terms of the rescaled parameters we get 


/i2 / 1 cosh(yL2^i(l - x)) \ ( 1 cosh(JLV2(l-x)) 

Vn,l[.x) = C„,1— - - +Cn,2— 

hi \L^hi L^sinhy^Zv / h2 yi 


h2fr 


Ly/JFi sinh J 


■ 9 {x) 


The parameters Cn^i and Cn ,2 from Eq. 39] are 


Uld ^n,l O'! 

Cn,l = fn{ - fbl) + ^-, - 

r /“2rf , , ^ , bn,l 02 

fn{ - + 62 ) + - 

e T h 2 ^ 


Cn,2 — 


with 


d = 


h2K2 L(3 


- y- - (616 + ^ 26 ) 


e = ai 6 + 026 H-f 


1 ^21 


Ki h2K-2 


and 


6 = 


coth coth \JL‘^h 2 


y/Jn 


6 = 




For a long DNA with L^hi ^ 1 we have 

h2 ( CnA 


VnM - 


: + 


L \hly/Jn ^ h2y/h^ 
, Cn,lZl2 , Cn,2h2 , fnh2 

W,i 0) = -+-+ —^ 

’ hi h2 P 

and from this we find for the sojourn times 


^n.2 


fj 


12 


3/3 




L 


Xn,l — 'f^n,l(0) T 

Un UnKl 


Lll2 ( Cn,l 


l^n yhly/hn h2y/h2j t'nKl V hi 


Xn,liO) 

Cn,2 \ ^ I-'ll2 / Cnp ^ Cn^2 
. . h2 


+ 


lufn 


(58) 


(59) 


(60) 

(61) 

(62) 




/? V 3 I)n XIJ 


(63) 


For ki 2 = 0 (/12 = 0) we recover riq = The number of switchings with yi = 00 are 


^12 — Tipki2 ~ LI‘^2 ( 


, + LHP 

\hly/Jn h2y/h2j 


hih2 

'W' 


(64) 
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Probability to bind to the target before detaching 

When a TF reaches the target in state 2 it eventually binds with probability P or it detaches 
with probability Q = 1 — P. To compute Q we consider a switching process with ^32 = 0 
(^32 = 0 ) to avoid rebinding to the DNA. The mean probability to detach before binding to the 
target when initially at the target site in state 2 is Q = /c 23 T' 2 , 2 ( 0 ) = hsf2^2(0)■ However, we 
cannot use f2,2(0) from our previous analysis because state 3 is missing and Vn,3{x) = 0 . More 
specifically, Eq. [ 23 ] is not valid when Vn,3{x) = I32 = 0 , and we cannot apply Eq. [ 37 |to obtain 
^22(0) as a function of ^21(0). We therefore recalculate ^21(0) and ^22(0) without Eg. (371 for 
I32 = 0 and the modified boundary condition f2,2(0) = The parameters a, / 3 , and fi2 are 
evaluated with I32 = 0 , e.g. /3 = I12I23 and a = li2 + hi + hs- By integrating Eq. [19] we find 


^2,1(0) — —^I2('r2,i — ^"2,2) 

^2,2(0) = 1 + ^21 (^2,1 — ^2,2) — 123X2,2 = 1 + (^21 + ^ 23 )(^ 2 ,l — 'f' 2 , 2 ) — 123X2,1 


(65) 


and from this we get 


, , ^21 + h3 , 1 

^^ 2 , 2 ( 0 ) + —;-^^2,i(0) - 1 


'^' 2 , 2 ( 0 ) + 1^2, 1 ( 0 ) “ 1 

tl2 



V. 


Ti(0) ^Q_ 


Kl 


23 


where we used 


= T2,2(0) = From 


( 66 ) 


(67) 


-r-X2,i{0) 

tl2 


/3 


— C2,lh2Zl — C2,2hlZ2 + C2,l -h C2,2- 

hi h2 

U2i(0) 

-C 2 ,lAi 26 - C2,2hl^2 H- ] -a - W'2(0) 


tl2 


= —C2,lh2^1 — C2,2hl^2 + 

■'l’2,l(0) = /23'72 ,i(0) — /2372,1 = 


^12 

^21 + h3 / I / 

- 1 -^2,1(0) + ^2,2(0) 

‘12 

h3 ^21 + hs 

Hi h2 


+ '‘’ 2 , 2 ( 0 ) — 1 


we get 


C 2 ,lh 2 il + C 2 , 2 hl ^2 = 1- —V 2 1 ( 0 ) 

Ki ’ 


( 68 ) 


Thus, we find the system of equations 


£ 2 , 11 ^ 2^1 + C 2 , 2 hi ^2 + —V 2 1 ( 0 ) 

Ki 


C2,lil + C2,2^2 + 


’^2,1 


( 0 ) 


Kl 


= 1 


9_ 

h3 


( 69 ) 
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By writing 


with 


we obtain 


C 2 ,l = + &1^^2,2(0) , C 2,2 = a2^^2,l(0) + &2^^2,2(0) 

h2 — hs 


hi ~ h'i 7 hi 

oi =-, 6i = 


hi - h2 


h2 - hi 


02 = 


h2 - hi 


b2 = 


h2 


hi - h2 


A 


4,1(0)] _ f 1 


4,2(o)J Vi^- 


with the matrix 


Aij — 


(aiiih2 + a2^2hi + ^ biCih2 + b2^2hi\ 

Y OlCl + 026 + ^ biCl + & 2?2 j 


The solution of Eq. [72] i 


IS 


^ 2 , 1 ( 0 ) 

^ 2 , 2 ( 0 ) 


1 


det(74) 

1 

det(^) 


A A ^ 
/±22 — - 


23 


— ^21 + ^11 


R 

I 23 


From this we obtain {P = 1 — Q) 

'^2,2(0) _ I23 t 


K2 = 


72,2(0) Q 
Q = I23 


= ^ 4 , 2 ( 0 ) = 


det(A) 




(70) 


(71) 


(72) 


(73) 


(74) 


(75) 


’ ^11 — K2det(^) 

For example, for k 23 = 0 or X 2 = 00 we have (5 = 0. The maximum is obtained for K 2 = 0 

^21 


Qmax — f 


23 


A 


11 


(76) 


Qmax depends on the switching rates and on ki. For example, Qmax = 1 is found for ki = 0 (no 
binding in state 1) or /21 = 0 (no switching to state 1). 











